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A geomètric framework for quantum statistical estimation is used to 
t***" ■ establish a series of higher order corrections to the Heisenberg uncertainty 

. relations associated with pairs of canonically conjugate variables. 

PACS Numbers : 02.40.K, 02.40.F, 02.50, 03.65 

\ An efíicient geomètric approach to parameter estimation, based upon Hilbert space ge- 

■ ometry, has been developed recently 0J. The idea is to regard a statistical model as being 
embedded into a real Hilbert space 7í, thereby allowing statistical estimation theory to be 
expressed neatly in terms of the geometry of TC. In particular, the concept of a statistical 
í> ! model can also be formulated in the case of a quantum mechanical state space, after a com- 
patible complex structure is introduced on the underlying real Hilbert space, thus leading 

O ' to a succinct characterisation of the quantum estimation problem. 
VO ' 
O 



In the present Letter, this geometrical framework is applied to establish a series of correc- 
. tions to the Heisenberg uncertainty relations. For example, in the case of the measurement 
of a 'position' parameter q, for which an unbiased estimator is given by an operator Q, 
with conjugate momentum operator P, where i[P, Q] — 1 (h — 1), we find an infinite set of 
corrections to the uncertainty lower bound, given in terms of higher order central moments 
of the distribution of P. If we write \x n = ((P — (P)) n ), these corrections take the form 
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The statistical model Ai in this case is the one-parameter family of states \ip q ) obtained 
by the action of the unitary operator e iqP on an initial state |^o)- The problem is to use 
measurements of the observable Q to estimate the value of q, provided that the actual state 
of the system lies somewhere on the manifold Ai. 

In order to pursue these ideas further, first we review a few of the bàsic concepts of 
statistical geometry. Consider a real Hilbert space 7í, with a symmetric inner product 
denoted g a b, and write £ a for a typical vector in TC. A random variable in Ti is represented 
by a symmetric operator X a b, with expectation = X a ^ a ^ b / '£ c £ c in the state £ a . The 

variance of X a f, is then 

Var ç [X] = ab \t , (2) 

where X a b = X a b — E^[X]g a b represents the deviation of X a b from its mean. In these 
expressions £ a is required to lie in the domain of X a b, which in the case of a bounded 
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operator is the whole of 7í. If £ a is normalised, we have g a b^ a ^ h = 1, whence £ a is a point 
on the unit sphere S in TC. By a 'statistical model', we mean a submanifold Ai of S, given 
parametrically by ^ a {9 1 ) where 9 l are local coordinates. The Riemannian mètric induced on 
Ai, known as the Fisher-Rao mètric, is given by [0,|J 

Qa = tgabdiÇ'd^ , (3) 

where <9j = d/89 1 (the factor of 4 is convent ional). For simplicity, in what follows we 
consider the case where Ai is one-dimensional, and for which the Fisher information is thus 
Ç = 4g ab £ a £ 6 , where the dot denotes d/ 89. 

Now we are able to pose the following statistical estimation problem. Given the mea- 
surement of an observable X abn we are interested in a one-parameter family £, a (9) of possible 
states characterising the distribution of the outcome. The parameter 9, which is under- 
stood to lie in some specified range, determines the unknown state of nature, and we wish 
to estimate 9 with the given data. For such an estimation problem a lower bound can be 
established for the variance of the estimate. Comparing the variance of our estimate to 
the lower bound, we can enquire to what extent the estimator is emcient. An elementary 
geometrical derivation of this bound is as follows. Given a curve £ a (#) in S, we say that a 
random variable is an unbiased estimator for a function r{9) if 

M = Ti0) • (4) 

along the given trajectory. Since r, regarded as a function of £ a , is homogeneous, we can 
extend its defmition to the whole of the domain of in 7í, and thus define the gradient 
V a T, where V a = d/d^ a . A short calculation shows that V a r = 1T ab ^ b / g cd ^ c ^ d , where 
T a è = T a b — rg a b, and as a consequence we deduce that the variance of the estimator T is 
given by 

Var ç [T] = ^ g ab V a rV b r (5) 

for any state on the given curve in S. This formula gives us a useful geometrical interpreta- 
tion for the variance of an estimator. 

The squared length of the gradient vector V a r can be expressed as the sum of the squares 
of its orthogonal components with respect to a suitable basis. To this end, we choose a set 
of vectors comprising the state £ a and its higher order derivatives. Letting £( n ) a denote the 
component of ^( n ) a = d n Ç a /d9 n orthogonal to both £ a and all of its derivatives of order less 

than n, we obtain a set of orthonormal vectors is given by {«f^V \J i^i^} {n = 0, 1, 2, • • •). 
It follows from the bàsic relation ([5]) that 

1 fí( n ) a V r) 2 
Var ç [T] > [ \ V °V • (6) 

Since each term in the sum on the right hand side of (0) is positive, we obtain a set of 
inequalities, which we refer to as generalised Bhattacharyya variance lower bounds 0. An 
essential difference here from the corresponding classical inequalities B is that the bounds in 



2 



(§) are not necessarily independent of the specific choice of estimator. In our approach to the 
quantum estimation problem below, however, we find a set of bounds that are systematically 
independent of the estimator T. 

We shall consider the case where M. is the submanifold of a quantum mechanical state 
space generated by the action of a one-parameter family of unitary transformations starting 
from some given initial state. In particular, we shall examine the case of state vector £ a (i) 
that satisfies the Schròdinger equation: 

t = r b H b c e ■ (7) 

Here, we identify 9 in the analysis above with the time parameter t, and the tensor J a b , satis- 
fying J a c J c b = —à a b , is the usual complex structure on 7i. The 'mean-adjusted' Hamiltonian 
H ab is symmetric and Hermitian, and describes the deviation of the Hamiltonian from its 
mean in the given state. The variance of the Hamiltonian is then given by AH 2 = g a b^ a ^ b - 

Note that in our approach the Hilbert space of quantum mechanics is treated as a real 
Hilbert space Tt endowed with a mètric g ab and a compatible complex structure J b . This is 
entirely equivalent to the conventional formulation of quantum mechanics based on complex 
Hilbert space U , but has the virtue of being more geometrically transparent in such a way, 
in particular, as to make the link with statistical geometry much clearer. 

Now, suppose that we have a large number of independent identical systems, each of 
which evolves from an initial state £ a (0) under the influence of the Hamiltonian H ab . The 
problem is to estimate how much time t has elapsed since the initial preparation. Let T ab 
be an unbiased estimator for the time parameter in this problem, so 

along the trajectory £ a (t). The existence of such an estimator is demonstrated, for exam- 
ple, in ||. As a consequence of (0) and @, it follows that 2T ab J b H^ a ^ d = 1 along the 
normalised trajectory £ a (í), which can be viewed as a 'weak' form of the canonical commu- 
tation relation. In other words, for an estimator, as opposed to an observable, we merely 
require an operator that satisfies the canonical commutation relation in expectation along 
the designated trajectory with reference to which the relevant measurement operation is 
being performed. 

We are interested in finding an explicit form of the variance lower bounds ([]) for quantum 
mechanical situations. When £ a satisfies the Schròdinger equation (0), we find that our series 
of orthogonal vectors is given by 

"ib t 

(-a ça { tbe \ ca S Çb l a 

Ç , ç > Ç — \s t,b)ç , ç "T^T - ^ ' ' ' 

Let {^/"} denote this series of vectors, so \l/g = £ a , ty® = £ a , and so on. We note that 
V tt í = 2T ab ^ b for normalised states in the domain of T ab , so we can rewrite @ in the form 

AT 2 > 1 E (TO 6 ) 2 , (9) 

" 4 V 9abK^ b n 
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Since can be expressed in terms of linear combinations of the derivatives of £ a , we find 
that for odd integers n, the following remarkable identity can be used to simplify equation 



2T ab ^ a e 



-irng ab e m)a c 
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1). Furthermore, a simple calculation shows that g a b£,^ a £,^ b = A t 2m, 



where m 

where ii% m = (H 2m ) denotes the 2m-th central moment of the Hamiltonian in the given 
state. Therefore, we find that for n = the contribution in (|D vanishes, and the n = 1 term 
gives (4/i 2 ) _1 , i.e., the Heisenberg lower bound. We also find |] that the n = 2 term depends 
upon the choice of the estimator T, while for n = 3 the result, which is independent of the 
choice of estimator, can be expressed in terms of the fourth cumulant of the distribution for 
the Hamiltonian H, as indicated in equation (1). We now proceed to deduce the general 
pattern of such higher order terms. 

In general, the even order contributions depend on the covariance of the estimator T 
with H 2n , while by virtue of (|Ï0|) the odd order terms are manifestly independent of the 
specific choice of T. Hence we only consider odd order terms. Let N n = g ab ^^t denote 
the denominator of the correction terms in (|9|). Then, an exercise shows that 



where D^n is defined by the determinant 
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therefore, D 2 = /i2, D 6 = HqH2 — H% and so on. We note, incidentally, that N 2 gives the 
curvature of the curve £ a (í) in S, and for k > 2 determines the higher order 'torsion' of 
£ a , in the sense of classical differential geometry. This can be seen by noticing that (O) is 
expressible in terms of the norm of the antisymmetric tensor obtained by skew-symmetrising 
the product of the following m + 1 vectors: $,j^~ 2 \ ■ ■ ■ , ^ ■ F° r example, we have 
D 10 = 6Çlakí 3) Çc\ ) Ç (5)[a Ç {3)b Ç {1)c] , and so on. Standard statistical identities § guarantee 
D211 > 0. Thus, in order to find the general pattern for the higher order corrections, we 
must understand the structure of the numerators 2\I/°T a ft£ 6 . Let us recali the way the basis 
vectors \P£ are formed: 









= e (3) 







That is, for each value of n, we subtract the components of from £( n ) a where k < n. 
Since we know 2^ a l f ah £ b = 1, we obtain 2^í afc £ 6 = /i^ 1 ( y u 4 - 3/í§) by use of (|Ï0]) and the 
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above expression for í^. Furthermore, from these two results, we can then find 2^/^T a bC, b , 
provided we know the coeficients of ^ (k = 1, 3) in ^5, and similarly for ^ and so on. 
In other words, once we find the pattern for these coefficients, we can express the resulting 
expression for higher order corrections in a recursive manner. 

We notice that the denominator terms in the coefficients are given by the expressions 
iVfc defined above. Hence we have reduced the problem to finding an expression for £^ a \l/fc a 
where k < n. If we define F n ^ = Ç( n ' a ty ka/^? 1 ^ kb, then after some àlgebra we deduce that: 



-1 



n-\-k _ 
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l^n+k ^n+k-2 ' ' ' A'n+l 
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In order to clarify the notation, we note that F n ^ and F n \ are given respectively by: 



n,3 



1 



-i\m+í 



A*n+1 
/l4 fi 2 



J^2 

with m = |(n — 1). We now have obtained all the relevant identities needed in order to find 
a recursive formula for the higher order corrections. Let us define the numerator in @ by 
U n = 2^/f l T a i ) ^ b . Then, from the following expression for 

K = - E F n, k n , (14) 

fe=l,3,5,- 

a recursive formula for U n can be obtained in the form: 

n-2 

U n = (-l) m n Mn _! - E F n,kU k , (15) 

fe=l,3,5,- 

with Ui — 1. Therefore, the uncertainty relation @ can be reexpressed as 

AT 2 AH 2 > - Y . (16) 

By input of the above relations (|Ï5|) and (|Ï6|), a suitable symbolic manipulation application 
can be made to produce expressions for the higher order terms purely in terms of the central 
moments fi 2 k of the distribution of the Hamiltonian H. 

One might wonder whether the higher order contributions in ( |T6"D vanish, e.g., for Gaus- 
sian distributions. However, it turns out that the expressions for U n in (|Ï5|) can be reex- 
pressed in terms of linear combinations of the cumulants n 2 k of the distribution for H. We 
recali that if 0(A) is the characteristic function (Fourier transform) of the distribution for 
H, then the cumulants defined by the relation 
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E«r^T- = ln< ^ (A) ' (17) 

and for Gaussian distributions all the cumulants higher than the second order vanish. 

On the other hand, for quantum mechanical applications, the distribution of H is not 
Gaussian, since H is bounded below. Hence one might wish to know typical orders of 
magnitude for these corrections. To obtain a crude estimate of the sort of numbers that 
might arise, we consider a system where the distribution of the Hamiltonian can be modelled 
by a gamma distribution, for which the density function is of the form 

p(H) = ^f aH H^ X , (18) 

where < H < oo and er, 7 > 0. In this case the resulting bound can be calculated, to 
second order, to be 

AT 2 AH 2 > 1 L + 18 \ 

- 4 ^ 3 7 2 + 47 7 + 42 ) ' y ! 

with A/í 2 = 7/a 2 . In particular, for 7 = 1, the distribution reduces to an exponential 
probability p(H) = a exp(-aH), and we find 

AT 2 AH 2 > _ (1 + 0.196 + 0.063 + •••) , (20) 

to third order, where the terms on the right are independent of a. These examples suggest 
an exponential decay for the higher order contributions, although in general this may not 
be the case. Extensive experimental data are available for the energy spectra of various 
systems, by use of which the associated higher order contributions can be evaluated. This 
would make for an interesting line of enquiry to pursue. One notices that in the case of 
an exponential distribution, the first two corrections already add up to a significant level of 
over 25% of the lower bound. 



We observe that all the corrections in (|I6) are given in terms of even order central 



moments of the distribution for H. This, however, does not imply that odd moments do not 
exist for quantum trajectòries. If the distribution for the conjugate variable is symmetric 
around its mean, then the odd moments vanish. However, as the situation stands, the system 
of orthogonal vectors does not form a complete set of basis vectors. In particular, on the 
unit sphere S, one can also consider the Cauchy-Riemann direction yf defined by r] a = J a b !; b - 
By virtue of the Schròdinger equation (0), i] a is automatically orthogonal, for example, to 
£ a and £ a . By use of this direction field and its derivatives, we can complete the basis. The 
inclusion of such terms in the expansion for the variance inequality leads to further refined 
bounds []7| , which are not generally negligible if the distribution is asymmetric. 

Our approach here has been to take a geomètric point of view in order to study the 
structure of the state space of quantum mechanics. The key to this approach is encapsulated 
in equation (|5|), where an essentially probabilistic, or statistical quantity, i.e., the variance, is 
expressed in terms of an essentially geomètric quantity, the length of a vector. Presumably, 
the higher order terms in the uncertainty relation ( |Ï6D can, in principle, be constructed by 
means of Standard Fourier analysis. In that respect it would be interesting to compare the 
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results we have obtained here to other approaches to generalised uncertainty relations ||. It 
should be emphasised, however, that the relations obtained therein take the form of 'strong' 
dispersion bounds, in the sense that they apply to all states in the Hilbert space (or at any 
rate a dense subset thereof); whereas our results are much more refined in the sense that 
they apply in a specific, physical measurement theoretic context for the determination of a 
parameter value in the case of a single trajectory of states (cf. ||). 

It is interesting to recali the fact that the higher order corrections can also be expressed by 
the cumulants of the distribution p(H). This result suggests the possibility of constructing 
a cumulant expansion for the uncertainty relations. The cumulants can also be associated 
with connected graphs with fixed numbers of 'legs', due to the logarithm in (|T7j), which can 
be compared with the use of Feynman diagrams for the free energy or ground state energy 
in quantum field theories. Therefore, it would be interesting to determine whether a neat 
graphic expansion for the uncertainty relations can be constructed. 

It should be pointed out that although the energy-time uncertainty is studied here as a 
matter of illustration, on account of its especially interesting features, analogous results hold 
for other conjugate variables, such as position and momentum, in which case the statistical 
manifold is given by the curve dÇ a /dq = J a b P b c ^ c for some initial state £ a (0). In this case 
both the momentum operator P and the position operator Q, which acts as an estimator for 
q, are self adjoint. The arguments outlined above then lead to the inequality ([1]). It should 
be born in mind that the geomètric formalism employed here, although highly effective, is 
not strictly essential, and the same results can be obtained by use of the familiar Dirac 
formalism, working in the Schròdinger representation. 
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